Determining the Saddle Point in Micromagnetic Models of Magnetization Switching 
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• A numerical model of single-domain nanoscale iron magnets fabricated using scanning- 

QQ ' microscope-assisted chemical vapor deposition is simulated using finite-temperature micromagnetics. 

A Projective-dynamics method is used to determine the magnetization at the saddle point as a func- 
' _i " tion of temperature. This magnetization is found to decrease linearly as the temperature is raised. 

^ ' An interesting problem in nonequilibrium statistical mechanics, with numerous applications in condensed-matter 

a physics and materials science, is how a system approaches the global minimum of the free energy. One simple, but 
common, situation is that of a local free-energy minimum which is separated from the global minimum by a free- 
energy barrier. In other words, any path between the metastable local minimum and the equilibrium global minimum 
involves an initial increase in the free energy. The most probable path between the minima involves the minimum free 
C/3 , energy increase since spontaneous increases in free energy are improbable. The maximum free energy along that most 
probable path corresponds to the saddle point. Often, the behavior of the nonequilibrium dynamics are dominated 
by properties near the minima and the saddle point. While the minima can be determined by examining histograms 
G of the state of the system, the saddle point is much harder to determine. In this paper we present a technologically 
important example of finding the saddle point using a projective-dynamics technique. 

The results described here are for a numerical model of single-domain nanoscale iron magnets that have been 
Q ' fabricated using chemical vapor deposition directed by a scanning-tunneling microscope These magnetic pillars 
O , have cross-sectional dimensions on the order of 10 nm and extend on the order of 100 nm perpendicular to the 
substrate. The magnetic pillars are modeled by a one-dimensional array of magnetization density vectors Ai{r) with 
fixed length Mg. Each vector precesses around a local field H{r) according to the Landau-Lifshitz-Gilbert equation 
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' where the electron gyromagnetic ratio 70 ~ 1.76 x 10 Hz/Oe. The phenomenological damping parameter a = 0.1 
, was chosen to give underdamped dynamics. Properties of bulk iron were assumed with the saturation magnetization 
• Ms — 1700 emu/cm'^ and the exchange length £x — 2.6 nm. The model pillar has a square cross section with area 
] and is SAi^ long. Details of the numerical approach, including the inclusion of thermal fluctuations in the local 

^ ■ field, are given in Ref. 0. 

I These nanopillars have a strong shape anisotropy, and the magnetization is most favorably oriented along the 
long axis of the magnet, which is taken to be the z direction. The nanomagnets are prepared in a metastable state 
via equilibration in an externally- applied magnetic field -\-Hoz, which is subsequently varied rapidly in magnitude 
Q . (without changing its orientation) to —Hqz. If the field Hq is less than the coercive field, here He ^ 1500 Oe a free- 
;V ■ energy barrier separates the positively-oriented magnetization from the equilibrium negative orientation. Previous 
simulations Q suggest that the saddle point in these pillars corresponds to the nucleation of a region of reversed 
magnetization at one of the ends of the pillar. The nucleated region grows until the entire magnetization is reversed. 

A useful technique for determining the magnetization corresponding to the free-energy maximum of the barrier, 
also called the saddle point, is the projective-dynamics technique ^ jj. The essence of this technique is projecting 
the original dynamics described in a high-dimensional phase space onto a probabilistic dynamic in a phase space of 
much lower dimension. A specific example is helpful. The dynamics of an Ising model can be projected onto the 
number of overturned spins, a measure of the global magnetization. The dynamics of the one-dimensional model are 
then described in terms of the probability that the number of overturned spins increases or decreases. Thinking in 
terms of the number of spins in the stable orientation, these probabilities are referred to as the probability of growing 
(Pgrow) and shrinking (Pshrink), respectively 1^. 

Projective dynamics for the model nanopillars is more complicated than that for the Ising model for two reasons. 
First, the magnetization is a continuous, three-dimensional variable. This is handled by projecting the ^-component of 
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the global magnetization into uniformly sized bins [[7| . Then Pgrow , the probability of the region of stable magnetization 
"growing," is the probability that during an integration step the magnetization moves into a bin corresponding to a 
smaller magnetization. Similarly, Pshrink is the probability that the magnetization moves into a bin corresponding to 
a larger magnetization. The second complication involves the persistent, subcritical regions of reduced magnetization 
that develop at each end of the pillar. To accomodate this, during the projective-dynamics analysis the pillar is divided 
into its top and bottom parts with the normalized, "global" magnetization of each half considered separately. 

The measured Pgrow and Pshrink at Hq—XQQQ Oe are shown in Fig. 1 for temperatures of 20 K, 50 K, and 100 K. 
For each temperature the results represent averaging over more than 10^ integration steps and on the order of 10'* 
switches. Consider the probabilistic dynamics at M2=0.85 for r=100 K. Here Pshrink>Pgrow and on average the 
magnetization moves to the right. For M^=0.98 the situation is reversed, and on average the magnetization moves to 
the left. The crossing point where Pshrink— Pgrow corresponds to a locally stable fixed point since, on average, here the 
magnetization moves towards the crossing. In this case the crossing is the metastable local free-energy minimum. For 
the crossing near Mz=0.74 the magnetization moves away from the point on average, and it corresponds to an unstable 
local maximum in the free energy. This is the saddle point. Measuring the value of Mz for which Pshrink= Pgrow is a 
convenient method for finding the saddle point. 

From Fig. 1 it is obvious that Mz at the saddle point is different for magnetization reversal at different temperatures. 
The value oi at the saddle point is presented vs temperature in Fig. 2. The point where Pshrink=Pgrow is estimated 
using lines fit via a least-squares method to each probability in the region near the crossing. The error bars are taken 
large enough to include nearly all for which the difference in the probabilities is less than the fluctuations in their 
estimates. 

There is a clear linear trend in at the saddle point with respect to the temperature. The solid line in Fig. 2 
is an unweighted least-squares fit to the data. This fit indicates that Mz{T — > 0)=0.860 at the saddle point. The 
dependence on temperature is quite strong, since Mz{100 K)=0.72. 

In summary, micromagnetic simulations at finite temperature have been used to investigate magnetization reversal 
in single-domain nanoscale magnets. The saddle point in the free energy was determined using the projective-dynamics 
technique, and the magnetization at the saddle point was found to depend linearly on the temperature for all the 
temperatures investigated, TfClOO K. 
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FIG. 1. Probability of "growing" and "shrinking" as functions of the magnetization of one half of the pillar for three 
temperatures: 100 K, 50 K, and 20 K. The left-most crossing shown for a given temperature indicates the saddle point, while 
the right-most crossing indicates the metastable local free-energy minimum. The crossings associated with the equilibrium 
global free-energy minimum lie near —1 and are not shown 
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FIG. 2. Magnetization along the pillar, Mz at the saddle point as a function of temperature (circle). The magnetization at 
the saddle point decreases with increasing temperature, and the solid line is a least-squaxes fit with intercept at 0.860. The 
error bars are estimated from the range of Mz where Pgrow andn Pshrink are nearly equal 
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